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CROSS-REFERENCE TO REALTED APPLICATION 

This application is a continuation-in-part of U.S. Application Serial No. 10/241,390, 
5 filed September 4, 2002. 

FIELD OF THE INVENTION 

The present invention is relates to the field of motor speed control. 

10 BACKGROUND 

Motor controllers govern the operation of a motor in order to direct activation, speed, 
torque, motor shaft position and other such characteristics of the motor during 
operation. Alternating current (AC) variable speed (motor) drives (VSD) may have 
encoderless feedback for control of the motor. Encoderless speed controls or 

15 estimators are based on the principle that the motor speed can be estimated from a 
measurement of the current or voltage waveform at a motor connection. Motor 
rotation induces other frequencies in the current waveform due to the physical 
construction (rotor bars or winding gaps) of the motor's rotor. A speed estimate can 
be generated by detecting these induced speed related frequencies and comparing 

20 them to known mathematical models relating the induced frequencies to the motor 
speed. Such a known mathematical model may be, for example, "The Nature of 
Induction Machines", P. Alger; Gordon and Breach, New York, 1965. 

SUMMARY 

25 In accordance with one aspect of the present invention there is provided a method of 
estimating a speed of a motor comprising: determining a correlation between a current 
wave sensed at the motor and frequency pairs from a set of weighted frequency pairs 
representing the current wave; fitting components of a motor control signal to a 
corresponding number of first orthogonal pairs in a set of weighted orthogonal pairs, 

30 the orthogonal pairs being orthogonal to the frequency pairs; fitting each of the 
frequency pairs as the orthogonal pair subsequent to the first orthogonal pairs; 
identifying the frequency pairs that provide a reduction in the mean squared error 
between the current wave and the set of weighted orthogonal pairs that satisfies a 
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criteria; and comparing desired frequencies from the identified frequency pairs with a 
motor speed harmonic model to determine an estimation of the speed. 

In accordance with another aspect of the present invention there is provided a method 
5 of estimating a speed of a motor using a current wave sensed at the motor, the method 
comprising: determining a correlation between a current wave sensed at the motor and 
frequency pairs from a set of weighted frequency pairs representing the current wave; 
fitting components of a motor control signal to a corresponding number of first 
orthogonal pairs in a set of weighted orthogonal pairs, the orthogonal pairs being 
10 orthogonal to the frequency pairs; fitting each of the frequency pairs as the orthogonal 
pair subsequent to the first orthogonal pairs; determining a reduction for each of the 
frequency pairs in a mean squared error between the current wave and the set of 
weighted orthogonal pairs according to a mean squared error between the current 
wave and the first orthogonal pairs; determining the reductions that satisfy a criteria 
15 and the corresponding frequency pairs; removing unwanted frequencies from the 

determined corresponding frequency pairs; and comparing the remaining frequencies 
from the determined corresponding frequency pairs with a motor speed harmonic 
model to determine an estimate for the speed. 

20 In accordance with a further aspect of the present invention there is provided a 

method of estimating a motor speed comprising: determining a correlation between a 
current wave sensed at the motor and frequency pairs from a set of weighted 
frequency pairs representing the current wave; fitting components of a motor control 
signal to a corresponding number of first orthogonal pairs in a set of weighted 

25 orthogonal pairs, the orthogonal pairs being orthogonal to the frequency pairs; 

comparing a subharmonic from the current wave with a harmonics speed model to 
identify regions in which to locate a corresponding harmonic, the subharmonic having 
a frequency less than a motor control signal; identifying a harmonics pair of 
frequencies in the regions having a separation from each other of no greater than a 

30 smallest harmonic of the motor control signal, wherein one of the frequencies is the 
harmonics pair of the corresponding harmonic; and comparing the harmonics pair 
with the harmonics speed model to determine an estimate for the speed. 
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In accordance with yet another aspect of the present invention there is provided a 
system for estimating a motor speed comprising: a correlation mechanism for 
determining a correlation between a current wave sensed at the motor and frequency 
pairs from a set of weighted frequency pairs representing the current wave; a fitting 

5 mechanism for fitting components of a motor control signal to at corresponding 

number of first orthogonal pairs in a set of weighted orthogonal pairs and fitting each 
of the frequency pairs as the orthogonal pair subsequent to the first orthogonal pairs, 
the orthogonal pairs being orthogonal to the frequency pairs; an mse reduction 
mechanism for determining a reduction in the mean squared error between the current 

10 wave and the set of weighted orthogonal pairs; a pair comparison mechanism for 
identifying the frequency pairs having a reduction that satisfies a criteria; a speed 
estimation mechanism for comparing desired frequencies from the identified 
frequency pairs with a motor speed harmonic model to determine an estimation of the 
speed; and a controller in communication with the correlation mechanism, the fitting 

15 mechanism, the mse reduction mechanism and the speed estimation mechanism for 
coordinating the process of estimating the motor speed. 

In accordance with an aspect of the present invention there is provided a system for 
estimating a motor speed comprising: a correlation mechanism for determining a 

20 correlation between a current wave sensed at the motor and frequency pairs from a set 
of weighted frequency pairs representing the current wave; a fitting mechanism for 
fitting components of a motor control signal to a corresponding number of first 
orthogonal pairs in a set of weighted orthogonal pairs, the orthogonal pairs being 
orthogonal to the frequency pairs; a region determination mechanism for comparing a 

25 subharmonic from the current wave with a harmonics speed model to identify two 
regions in which to locate a corresponding harmonic, the subharmonic having a 
frequency less than a motor control signal; a corresponding frequencies mechanism 
for identifying a harmonics pair of frequencies in the two regions having a separation 
from each other no greater than a smallest harmonic of the motor control signal, 

30 wherein one f the frequencies in the harmonics pair is the corresponding harmonic; a 
speed estimation mechanism for comparing desired frequencies from the identified 
harmonics pairs with a harmonics speed model to determine an estimation of the 
speed; and a controller in communication with the correlation mechanism, the fitting 
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mechanism, the mse reduction mechanism and the speed estimation mechanism for 
coordinating the process of estimating the motor speed. 

In accordance with an aspect of the present invention there is provided a computer 
readable medium having stored thereon computer-executable instructions for 
estimating a speed of a motor comprising: determining a correlation between a current 
wave sensed at the motor and frequency pairs from a set of weighted frequency pairs 
representing the current wave; fitting components of a motor control signal to a 
corresponding number of first orthogonal pairs in a set of weighted orthogonal pairs, 
the orthogonal pairs being orthogonal to the frequency pairs; fitting each of the 
frequency pairs as the orthogonal pair subsequent to the first orthogonal pairs; 
identifying the frequency that provide a reduction in the mean squared error between 
the current wave and the set of weighted orthogonal pairs that satisfies a criteria; and 
comparing desired frequencies from the identified frequency pairs with a motor speed 
harmonic model to determine an estimation of the speed. 

In accordance with an aspect of the present invention there is provided a computer 
readable medium having stored thereon computer-executable instructions for 
estimating a motor speed comprising: determining a correlation between a current 
wave sensed a the motor and frequency pairs from a set of weighted frequency pairs 
representing the current wave; fitting components of a motor control signal to a 
corresponding number of first orthogonal pairs in a set of weighted orthogonal pairs, 
the orthogonal pairs being orthogonal to the frequency pairs; comparing a 
subharmonic from the current wave with a harmonics speed model to identify regions 
in which to locate a corresponding harmonic, the subharmonic having a frequency 
less than a motor control signal; identifying a pair of frequencies in the regions having 
a separation from each other no greater than a smallest harmonic of the motor control 
signal, wherein one of the frequencies in the pair is the corresponding harmonic; and 
comparing the pair of frequencies with the motor speed harmonic model to determine 
an estimate for the speed. 

In accordance with an aspect of the present invention there is provided a computer 
readable medium having stored thereon computer-executable instructions for 



estimating a speed of a motor using a current wave sensed at the motor, the method 
comprising: determining a correlation between a current wave sensed at the motor and 
frequency pairs from a set of weighted frequency pairs representing the current wave; 
fitting components of a motor control signal to a corresponding number of first 

5 orthogonal pairs in a set of weighted orthogonal pairs, the orthogonal pairs being 

orthogonal to the frequency pairs; fitting each of the frequency pairs as the orthogonal 
pair subsequent to the first orthogonal pairs; determining a reduction for each of the 
frequency pairs in a mean squared error between the current wave and the set of 
weighted orthogonal pairs according to a mean squared error between the current 

l o wave and the first orthogonal pairs; determining the reductions that satisfy a criteria 
and the corresponding frequency pairs; removing unwanted frequencies from the 
determined corresponding frequency pairs; and comparing the remaining frequencies 
from the determined corresponding frequency pairs with a motor speed harmonic 
model to determine an estimate for the speed. 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention will be described in conjunction with the drawings in which: 
Fig. 1 illustrates a motor control system according to an embodiment of the present 
invention; 

20 Figs. 2A-B illustrate a method of estimating the speed of a motor according to a first 
embodiment of the present invention; 

Fig. 3 illustrates a speed estimator of the motor control system according to an 
embodiment of the present invention; 

Fig. 4 illustrates the frequency spectrum in which the speed estimator detects speed 
25 related information; 

Figs. 5A-B illustrate a method of estimating the speed of a motor according to a 
second embodiment of the present invention; and 

Fig. 6 illustrates a speed estimator of the motor control system according to the ) 
second embodiment. 

30 

DETAILED DESCRIPTION 

Figure 1 illustrates a motor control system 100 according to an embodiment of the 
present invention. A vector controller 102 in the motor control system 100 drives an 
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AC induction motor 104 via an inverter 1 10 and motor drive connections 112. The 
motor 104 is driven by a known frequency, the fundamental frequency, which is 
produced by the vector controller 102. The vector controller 102 receives information 
from a rotor flux reference 106, a speed reference 108 and an estimated speed of the 
rotor of the motor 104 from a speed estimator 114. Based on this information the 
vector controller 102 provides control signals 124 to the inverter 110 for control of the 
flow of power from a power source 1 16 to the motor 104. As the inverter 110 
receives the control signals 124 from the vector controller 102, power from the power 
source 1 16 to the motor 104 may be regulated thereby. 

The AC induction motor 104 includes a stator and the rotor (not shown) that rotates in 
response to the power provided to the motor 1 04. Current induced in the rotor as a 
result of a magnetic field produced by power supplied to the stator includes 
frequencies representative of the rate of rotation of the rotor and frequencies produced 
during rotation due to the physical construction of the rotor (e.g. rotor bars or winding 
gaps, etc.). This current is measured by sensors 122 at the motor drive connections 
1 12 and forms the current wave from which the speed estimator 114 derives an 
estimate of the motor speed. 

The rotor flux reference 1 06 is a measure of the rotor flux-linkage and provides a 
reference frame for measurement of torque in the motor 104. This reference frame 
can be used to control the torque by regulating the rotor flux-linkage and stator 
current. 

The speed reference 108 is an indication of the desired operating speed of the rotor. 
The vector controller 102 uses the speed reference 108 as a baseline against which the 
motor speed obtained from the speed estimator 1 14 can be compared. The control 
signals 124 produced by the vector controller 102 are generated to control the motor 
1 04 to align the motor speed with the speed reference 1 08. That is, the fundamental 
frequency, also called the carrier drive frequency (CDF), or the amplitude of the CDF 
produced by the vector controller 102 is generated to adjust the speed of the motor 
1 04 so that it is the same as the speed reference 1 08. 



A voltage feedback may be provided from the motor drive connections 1 12 to the 
vector controller 102. 

The vector controller 1 02 provides the speed estimator 1 14 with the CDF 1 1 8. The 
5 vector controller 1 02 can be based on any of a number of known vector control 
methods, such as the one set forth in "Vector Control of AC Machines," P. Vas; 
Oxford Science Publications, 1990. For a synchronous motor a motor speed 
controller similar to the vector controller 102 may be used. 

10 The speed estimator 1 1 4 determines an estimate of the motor speed based on samples 
of the current wave sensed by the sensors 122 at the motor drive connections 1 12 and 
the CDF 1 1 8. For low speed operations, the speed estimator 1 14 alternatively 
receives a selected carrier drive frequency (SCDF) 120 from the vector controller 102. 

1 5 The CDF 1 1 8 provided to the speed estimator 1 1 4 is the carrier drive frequency 
produced by the vector controller 1 04 at the time when the current samples where 
obtained from the motor drive connections 112. In order to provide a constant CDF 
during current wave obtainment, the speed estimator 114 may inform the vector 
controller 102 that the current at the motor drive connections 1 12 is being sampled. 

20 In response, the vector controller 1 02 may produce a constant CDF for a 

predetermined amount of time during which the speed estimator 1 14 obtains current 
wave samples from the sensors 122. 

When the CDF values are small and the motor 1 04 is operating at only a small 
25 fraction of the base speed then the CDF 1 1 8 and the SCDF 1 20 may form the control 
signal 1 24. The SCDF 1 20 is of a higher frequency but lower amplitude than the 
CDF 118. The CDF 1 1 8 and the SCDF 120 are added together to produce speed 
related harmonics in the current wave from the sensors 122 to be detected by the 
speed estimator 114. The SCDF 120 is selected such that the speed related harmonics 
30 have a sufficient separation for detection. Although adding the SCDF 1 20 decreases 
the motor's efficiency due to waste of some power, the SCDF 1 20 enables the speed 
estimator 1 14 to estimate low motor speeds where the frequencies induced by the 
CDF 1 18 are closely spaced together (e.g. 1 Hz) making it hard to distinguish them. 
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The SCDF 120 may be a preset frequency used regardless of the speed (e.g. 30 Hz or 
60 Hz for all occasions that the SCDF 120 is used). 

The speed estimator 1 1 4 uses the CDF 1 1 8/SCDF 1 20 and the sampled current wave 
in a Real Time Fast Orthogonal Search (RTFOS) process to estimate harmonics 
representative of the motor speed in order to produce the motor speed estimate. 

The inverter 1 10 receives the control signals 124 from the vector controller 102 to 
control flow of power from the power source 1 16 to the motor 104. When the speed 
received from the speed estimator 1 14 does not meet certain minimum criteria, the 
vector controller 102 enhances the CDF 1 1 8 by adding the SCDF 120 thereto. The 
combined CDF 1 1 8 and SCDF 1 20 form the control signal 1 24 that is provided to the 
inverter 1 1 0. The inverter 1 1 0 may be formed of an array of field effect transistors 
that couple power from the power source 1 1 6 to the motor 1 04. 

Figs. 2A-B illustrate a method 200 of estimating motor speed by the speed estimator 
1 14 according to an embodiment of the present invention. The method 200 is based 
on the Fast Orthogonal Search (FOS) method. FOS uses a set of non-orthogonal 
functions (e.g. a portion of a sine and a portion of a cosine wave) and fits them to a 
sampled signal (e.g. the current wave). FOS searches for frequencies with a period 
longer than a specific sampling window as well as those with a fractional number of 
periods in the sampling window. 

When the motor 104 is in operation current is induced in the motor drive connections 
1 12 from the rotor of the motor 1 04. If the CDF is small then the vector controller 
102 may add the SCDF to the CDF, thus the SCDF may also be received in step 204. 
The sensors 122 connected to the motor drive connections 1 12 detect and measure 
this current wave. The vector controller 1 02 is informed in step 202 that the current 
wave from the motor drive connections 1 1 2 is going to be measured so that the vector 
controller 102 can produce a constant CDF. 

The measured current wave is received from the sensors 1 22 and the CDF is received 
from the vector controller 102 in step 204. The measured current wave is converted 
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to a digital sampled current wave in step 206. If the sensors 122 detecting the current 
wave provide a sampled current wave that is already digital then step 206 need not be 
performed. A predetermined number of samples N of the current wave are measured 
in step 208. Each of the N samples may be taken such that each sample is separated 

5 by a fixed amount of time from consecutive samples. The amount of time spent 
collecting the samples may depend on the number of speed estimates per second as 
limited by the ability of the speed estimator 1 14 to estimate the speed and store the 
estimated speed or measured samples. Accuracy of the vector controller 102 (i.e. 
speed regulator) improves with an increase in the number of speed estimates per 

10 second that are measured. If 6 to 1 0 speed estimates per second are produced by the 
speed estimator 1 14 then control of the motor 104 operating below 2% of the base 
rated motor speed will generally be stable. 

The current wave samples can be modeled as a weighted sum of candidate frequency 
15 pairs by: 

M 

y(n) = ^a mPm (n) + <;(n) (D 



m=0 



( 27if n\ ( 1rf m n\ 

where p 2m (n) = sin| -^f- and p 2m+] (n) = cosl -y- 



(2) 



where y(n) is the sampled current wave which is a weighted (a m ) combination of the 
candidate functions p m (n) and the error ^(n) between y(n) and the functional 
20 expansion (sum of the weighted candidate functions). The candidate function pairs 
p 2m (n) and p 2m +i(n) representing the candidate function are based on the sampling 
frequency F s and the set of candidate frequencies f m . The candidate frequency pairs 
are chosen to represent frequencies between 0Hz and the CDF. 

25 The sampled current wave is correlated to a set of candidate frequency pairs in step 
210. The candidate frequency pairs are chosen to represent sine and cosine function 
pairs (e.g. the candidate functions) having particular frequencies of interest and 
remain the same every time the method 200 is executed. The sampled current wave 
can be represented by a weighted combination of candidate functions and an error. 

30 The candidate frequency pairs represent the candidate function from which the 

sampled current wave can be derived. A subset of candidate frequency pairs that is 
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noncontiguous can be used as a means of reducing the number of calculations 
performed by the speed estimator 1 14. 



Correlation of the current wave with the set of candidate frequency pairs may be 
performed via known correlation techniques or by using zero-filled, adjusted Fast 
Fourier Transform (FFT). In the case of the later correlation calculation, the set of 
candidate frequency pairs are chosen to correspond to frequency bins of the zero- 
filled FFT. Zeros are added to the sampled signal to increase their length, such as by 
adding 1024 zeros to the end of 1024 sample values to increase the set length to 2048. 
Zero-filling increases the number of candidate frequencies calculated for correlation 
by the FFT. If it is known that that the FFT is zero-filled, some of the multiplication 
and addition operations in the FFT algorithm can be removed. 

The candidate frequency pairs are determined to correspond with the FFT bins and to 
have a higher frequency resolution than the FFT. The candidate frequency pairs are 
separated from each other by approximately four times the FFT bin size. 

Correlation of the sampled current wave and the set of candidate frequency pairs may 
be performed using a modified Fast Fourier Transform (FFT) given by : 



N-\ 



2nk 



Y(k) = j^y(n)e~ M ( 3 ) 

»=0 

where k is the frequency bin index, n is the time index and M is the number of 
samples in the time series after it has been zero-filled. This complex equation (3) can 
be expressed as sinusoids and cosinusoids resulting in: 

y ( ,)=X ,(")co(— J-,Z *»"\-JT) . (4) 



The sampled current wave FFT is zero-filled to M points to increase the number of 
frequency bins calculated. Since the candidate frequencies should coincide with FFT 
bin frequencies and the resolution of the FFT is R FFT =Fs/M then 

f = (5) 

m M 

where 0<k<M/2. 
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A suitable type of FFT is the real valued FFT with Decimation in Frequency (DIP). 
The DIF FFT processes the input in the order collected. Since the bottom half or 
three-quarters of the input is zero, the multiplication and addition operations required 
by these zero elements need not be performed, thus reducing the number of 
computations performed by the method 200. 

While the sampled current wave y(n) can be represented as a weighted combination of 
candidate functions p m (n), the sampled current wave y(n) can also be presented by a 
combination of weighted orthogonal functions w m (n) and some error e(n) as follows: 

v(H) = i> m w m (n) + *(«) ( 6 ) 



m=0 



where g m are the weights for the orthogonal functions. 

When the weights g m relating the orthogonal functions w m (n) to the sampled current 
wave y(n) are determined they can be used to determine the weights a m relating the 
non-orthogonal candidate functions p m (n) to y(n). The orthogonal functions are 
derived from the candidate functions and may be represented by orthogonal pairs 
similar to the candidate frequency pairs. 

Such derivation of the orthogonal functions may be performed using a modification of 
the Gram-Schmidt orthogonalization algorithm. The Gram-Schmidt algorithm is 
described in "Fast Orthogonal Search for Array Processing and Spectrum 
Estimation," K. Adeney and M. Korenberg; IEE. Proc. Visual Image Signal 
Processing, Vol. 141, No. 1 Feb. 1994, herein incorporated by reference. 

If the CDF is not less than 1Hz or the motor speed is greater than 2% base rated speed 
as determined in step 212 then the CDF is fit as the first candidate frequency pair in 
the orthogonal pairs set in step 214. 

The process of fitting the candidate frequency pairs involves deriving a function 
orthogonal to all previously fitted functions based on the candidate frequency pair. 
As the orthogonal functions are derived, the weights of these orthogonal functions are 
determined as part of the fitting process. 
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An orthogonal function can be determined using the corresponding candidate function 
minus the weighted sum of previous orthogonal functions, as given by: 

m-1 

w m (n) = ^ m («)-Z flf -^(") (7) 

where a is an m x m matrix containing for Gram-Schmidt weights used to derive the 
orthogonal functions from the previous orthdgonal fiinctions Only half of the matrix 
a is used thus only the non-zero elements are stored thereby reducing the memory 
used to store a. 

The first orthogonal function is set equal to the first candidate pair (the CDF). 

w 0 (n) = i? o (n) ( 8 ) 

A time correlation function is defined as: 

D(m,r) = p m (n)w r {n) (9) 
and it can be shown that : 



D(m,m) = w 2 m (n) 0°) 
The weight a™, relating w m to w r can be computed by : 



P„,("M(") = £j^hll (11) 



w r 2 (n) W> 



Since determining the orthogonal functions w n (n) is not required, as the weights g m 
are used to determined the weights a m of the non-orthogonal functions, D(m,r) can be 
determined as follows: 

r-1 

£>(/», r) = p m (n)p r (n)-Y, a rMm,i) (12) 



The weights g m of the orthogonal series is: 

c ( m ) (13) 



Sm D(m,m) 



where 
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C{m) = y{n)p m {ri)-Y, a m rC{r) 

r=0 



(14) 



It can be shown that the weights relating the candidate functions to the sampled signal 
are given by: 

i-m 

where 

v. = -f> mr v r i=m+l, m+2, ...,M (16) 

r=m 

If the CDF is less than 1 Hz and the motor speed is less than 2% rated base speed as 
determined in step 21 2 then the CDF is fit as the first candidate frequency pair and the 
SCDF is fit as the second candidate frequency pair in the orthogonal set in step 216. 

After the first candidate frequency pair(s) has been fit in the orthogonal pairs in either 
step 214 or step 216, the mean squared error (MSE) between the orthogonal functions 
set in step 214 or step 216 and the sampled current wave is determined in step 218. 
The method 200 attempts to minimize the mean squared error between the orthogonal 
functions and the sampled current wave. 

The mean squared error is given by: 

M 

* 2 («) = v 2 («)-!>>>) ( g ) 

m=0 

The term g 2 m w„ 2 ,(n) is the mean squared error reduction caused by adding the mth 
candidate frequency pair to the orthogonal set. Since each frequency component of 
the sampled current wave is composed of a sine and a cosine component, the mean 
squared error reduction for both the sine and cosine candidate of a particular 
frequency is determined and combined to form the overall mean squared error 
reduction for that frequency f m . 

Each candidate frequency pair in the set is fitted as the second frequency pair (or third 
if the SCDF is used) in the set orthogonal pairs in step 220 in a manner similar to the 
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fittings in steps 216 and 214. Note that since the frequency corresponding to the CDF 
(or CDF and SCDF) has already been fit it will not be in the remaining candidate 
frequency pairs. The means squared error reduction for each additional pair in the 
orthogonal sets resulting from step 220 is determined in step 222. 

After all of the candidate frequency pairs have been tried as the second pair (or third 
pair if the SCDF was used) in the orthogonal set and the resulting mean squared error 
reduction determined, a predetermined number of candidate frequency pairs providing 
the highest reduction in the mean squared error are selected in step 224. 
Alternatively, all candidate pairs providing a reduction within a predetermined range 
of the highest mean squared error reduction may be selected. Selecting more than one 
candidate frequency pair at a time reduces the number of computations performed by 
the method 200. 

The selected candidate frequency pairs are examined in step 226 to remove any pairs 
that are CDF harmonics. This may be done by looking for frequencies separated by 
specific amounts. The remaining selected candidate frequency pairs are inserted into 
eh orthogonal pairs in step 228 

Steps 220 to 228 are repeated until a specified criteria is met as determined in step 
230. The specified criteria may include fitting a determined number of candidate 
pairs, fitting enough candidate pairs to reduce the mean squared error to a certain 
threshold, continue to fit candidate pairs until there is no reduction in the mean 
squared error. 

The remaining pairs are assumed to be motor speed harmonics. A speed estimate can 
be generated from the remaining selected candidate frequency pairs by comparing 
them to a known mathematical model relating the induced frequencies to the motor 
speed in step 232. An exemplary mathematical model suitable for use is described in 
"The Nature of Induction Machines," P. Alger; Gordon and Breach, New Work, 1965. 

Fig. 3 is a system diagram showing the speed estimator 1 14. The speed estimator 1 14 
includes a controller 304 that interconnects a sampling mechanism 306, a correlation 
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mechanism 308, a fitting mechanism 310, an MSE mechanism 314, a vector 
controller interface 324, a CDF analysis mechanism 312, an MSE reduction 
mechanism 316, a pair comparison mechanism 318, a harmonics recognition 
mechanism 320 and a speed estimation mechanism 322. 

The controller 304 coordinates processing between these modules including 
facilitation of data movement and flow control. 

The vector controller interface 324 interfaces with the vector controller 102 to send an 
indication that the speed estimator 1 14 is going to measure the current wave in the 
motor drive connections 112 thus allowing the vector controller 102 to send a 
constant CDF for a predetermined period of time. The CDF (or SCDF and CDF) that 
is produced by the vector controller 102 during this time is received at the speed 
estimator 1 1 4 by the vector controller interface 324. 

A sensor interface 300 is connected to the sampling mechanism 306 through an A/D 
converter 302. The sensor interface 300 is in communication with the sensors 122 at 
the motor drive connections 1 12 to receive the current wave measured therefrom. The 
measured current wave received at the sensor interface 300 is provided to the A/D 
converter 302 for conversion from an analog signal to a digital signal, thus forming 
the sampled current wave. The key information in the current wave is about 60db 
smaller than the CDF. 

The A/D converter 302 may be any known A/D converter capable of handling such 
data. If the current wave received from the sensors 122 is already a digital signal then 
the current wave need not be processed by the A/D converter 302. 

The sampling mechanism 306 receives the sampled current wave and extracts N 
samples therefrom to produce the sampled current wave that will be used to estimate 
the speed of the motor 104. 

The correlation mechanism 308 receives the sampled current wave consisting of the N 
samples and determines a correlation between the sampled current wave and a set of 
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candidate frequency pairs that are produced by the controller 304. The sampled 
current wave may be represented as a weighted sum of a subset of candidate 
frequency pairs. 

The CDF analysis mechanism 312 determines whether the CDF from the vector 
controller 1 02 is less than 1Hz and the motor speed is less than 2% rated base speed, 
indicating low speed operation of the motor 104. If the above holds true then the 
SCDF is used in addition to the CDF. That is, the CDF analysis mechanism 312 
classifies a previous motor speed or the speed reference 1 08 as being either low or 
normal. 

The fitting mechanism 310 fits the candidate frequency pairs in an orthogonal set, 
which may also be used to represent the sampled wave current. The first pair that is 
fit in the orthogonal set is the CDF. If the SCDF was used then the SCDF is fit as the 
second pair in the orthogonal set. The fitting mechanism 310 comprises a frequency 
weights mechanism 326 and an orthogonal weights mechanism 328. The orthogonal 
weights mechanism 328 determines the weights for the orthogonal set based on the 
values of the frequency pairs. The frequency weights mechanism 326 determines 
weights for the candidate frequency pairs based on the orthogonal weights. 

The MSE mechanism 314 determines the mean squared error between the orthogonal 
set and the sampled current wave. The MSE reduction mechanism 316 determines the 
mean squared error reduction when additional candidate frequency pairs are fit into 
the orthogonal set. The pair comparison mechanism 318 determines those pairs that 
provide the highest mean squared error reduction between the sampled current wave 
and the orthogonal set when they are added to the orthogonal set. The pair 
comparison mechanism 318 may select a predetermined number of candidate 
frequency pairs that provide the highest mean squared error reduction or may select 
all candidate pairs providing a mean squared error reduction within a determined 
range of the highest value. 

The harmonics recognition mechanism 320 determines those candidate frequency 
pairs from the selected candidate frequency pairs determined by the pair comparison 



16 



mechanism 318 that are harmonics of the CDF and removes them from the selected 
candidate frequency pairs. 



The speed estimation mechanism 322 estimates the speed of the motor 1 04 based on 
the selected candidate frequency pairs remaining after the CDF harmonics have been 
removed. The speed estimation is based on a comparison of the remaining selected 
candidate frequency pairs with a known mathematical model relating the induced 
frequencies in the motor drive connections 112 with the motor speed. 

Fig. 4 depicts a frequency spectrum graph of the measured current wave, filtered and 
matched for a speed estimate by the speed estimator 114. The sampling sensitivity of 
the A/D converter 302 is shown 406 with reference to the vertical (amplitude) axis. In 
this example, the CDF is 60Hz and drives the motor 104 at the rated base speed 401 
with resulting unwanted odd harmonics 404 (even and other harmonics of similar 
amplitude are not shown for clarity). A speed-induced sub-harmonic frequency 402 is 
shown at 25Hz that indicated a rough estimate of motor speed as: 

25Hz/CDF (or CDF and SCDF) (60Hz) x rated base speed = 750 rpm 
where the example rated base speed is 1 800RPM. Other unrelated noise frequencies 
are shown to either side 403 of the sub-harmonic. Examples of two speed-induced 
harmonics are shown at 120Hz and 1 80Hz 405 separated by 2 times the CDF. Note 
that the second frequency, 1 80Hz, coincides with a harmonic of the CDF; 
consequently, the frequency may not be used in the final speed estimate 
determination. 

Figs. 5A-B illustrate a speed estimation method 500 according to a second 
embodiment of the present invention. Steps 502 to 518 of the method 500 
corresponding with step 202 to 21 8 of the method 200 of the first embodiment and 
will not be discussed in detail. 

After the mean squared error between the orthogonal set containing the CDF (or CDF 
and SCDF) and the sampled current wave has been determined, a search for possible 
speed-induced sub-harmonic frequencies is performed between the CDF (or SCDF) 
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and OHz in step 520. Sub-harmonic frequencies are those frequencies that are speed- 
induced frequencies harmonics but have a lower frequency than the CDF (or SCDF). 

Each speed-induced sub-harmonic determined in step 520 is compared to a model 
relating the induced frequencies with the estimated speed. The purpose of this 
comparison is to locate two separate frequency regions where corresponding speed 
induced harmonics would have to be in order to make a pattern match with the speed 
model in step 522. These frequency regions correspond with the FFT bins used in the 
correlation in step 510. These regions may be equivalent to 1 FFT bin higher and 1 
FFT bin lower than the identified sub-harmonics. 

Any frequency harmonics of the sub-harmonics in each identified frequency region 
are then located in step 524. During a search of each region various techniques may 
be used to filter and/or discriminate against other noise carriers using learned methods 
to avoid the time wasted detecting and processing erroneous signals (e.g. an inversion 
of the induced harmonics). 

The speed-induced frequency harmonics found are examined to determine if the 
harmonic frequencies correspond to a second harmonic frequency whose frequency is 
two times the CDF away in step 526 or two harmonic frequencies that are separated 
by the smallest harmonic of the control signal (e.g. CDF or CDF+ SCDF). 

In step 528 it is determined whether there were any frequency pairs determined in step 
526. If there was at least one frequency pair was determined then the frequency pair 
that has the greatest separation from the location of known harmonics (e.g. CDF 
harmonics) is selected in step 530 to avoid unwanted motor rather than speed related 
frequencies. The selected frequency pair is compared in step 532 to the known 
harmonic speed model to estimate the motor speed. If there is more than one valid 
speed estimate, the estimates are averaged to result in a final value to drive the vector 
controller 104 in step 534. 

If there were no frequency pairs determined in step 528 then it is determined in step 
536 if other sub-harmonics were found. If other sub-harmonics were found then steps 
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522 to 528 are repeated. If no other sub-harmonics were found then no new speed 
estimate is produced in step 538. 



Matching the two frequencies as in step 526 is not the only means of separating 
speed-induced harmonics from unwanted harmonics. Frequency pairs that coincide 
with known harmonics may be discarded. These steps may be replaced with more 
sophisticated pattern matching algorithms and a speed tracking routine. 

Fig 6 is a system diagram showing the speed estimator 114 according to the second 
embodiment. The speed estimator 1 14 includes a controller 604 that interconnects a 
sampling mechanism 606, a correlation mechanism 608, a fitting mechanism 610, an 
MSE mechanism 614, a vector controller interface 624, a CDF analysis mechanism 
612, a sub-harmonics mechanism 615, a region determination mechanism 618, a 
harmonics identification mechanism 61 0 a corresponding frequencies mechanism 626 
and a speed estimation mechanism 622. 

The controller 604 coordinates processing between these modules including 
facilitation of data movement and flow control. 

The vector controller interface 624 interfaces with the vector controller 102 to send an 
indication that the speed estimator 1 14 is going to measure the current wave in the 
motor drive connections 112 thus allowing the vector controller 102 to send a 
constant CDF for a predetermined period of time. The CDF (or SCDF and CDF) that 
is produced by the vector controller 1 02 during this time is received at the speed 
estimator 1 14 by the vector controller interface 624. 

A sensor interface 600 is connected to the sampling mechanism 606 through an A/D 
converter 602. The sensor interface 600 is in communication with the sensors 122 at 
the motor drive connections 1 12 to receive the current wave measured therefrom. The 
measured current wave received at the sensor interface 300 is provided to the A/D 
converter 602 for conversion from an analog signal to a digital signal, thus forming 
the sampled current wave. The key information in the current wave is about 60db 
smaller than the CDF. 
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The A/D converter 602 may be any known A/D converter capable of handling such 
data. If the current wave received from the sensors 122 is already a digital signal then 
the current wave need not be processed by the A/D converter 6302. 

The sampling mechanism 606 receives the sampled current wave and extracts N 
samples therefrom to produce the sampled current wave that will be used to estimate 
the speed of the motor 104. 

The correlation mechanism 608 receives the sampled current wave consisting of the N 
samples and determines a correlation between the sampled current wave and a set of 
candidate frequency pairs that are produced by the controller 604. The sampled 
current wave may be represented as a weighted sum of a subset of candidate 
frequency pairs. 

The CDF analysis mechanism 612 determines whether the CDF from the vector 
controller 102 is less than 1Hz and the motor speed is less than 2% rated base speed, 
indicating low speed operation of the motor 1 04. If the above holds true then the 
SCDF is used in addition to the CDF. That is, the CDF analysis mechanism 312 
classifies a previous motor speed or the speed reference 108 as being either low or 
normal. 

The fitting mechanism 610 fits the candidate frequency pairs in an orthogonal set, 
which may also be used to represent the sampled wave current. The first pair that is 
fit in the orthogonal set is the CDF. If the SCDF was used then the SCDF is fit as the 
second pair in the orthogonal set. The fitting mechanism 61 0 comprises a frequency 
weights mechanism 630 and an orthogonal weights mechanism 628. The orthogonal 
weights mechanism 628 determines the weights for the orthogonal set based on the 
values of the frequency pairs. The frequency weights mechanism 630 determines 
weights for the candidate frequency pairs based on the orthogonal weights. 

The MSE mechanism 614 determines the mean squared error between the orthogonal 
set and the sampled current wave. 
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The sub-harmonics mechanism 616 locations all speed-induced sub-harmonics in the 
current wave between 0Hz and the CDF (or SCDF and CDF). 

The region determination mechanism 618 compares one of the determined sub- 
harmonics with a harmonic speed model to identify frequency search regions. 

The harmonics identification mechanism 620 locates speed-induced harmonics of the 
sub-harmonics in the identified frequency search regions in the current wave. 

The corresponding frequencies mechanism 626 determines frequencies in one of the 
regions determined by the region determination mechanism 61 8 that have a 
corresponding frequency in the second region where there is a spread of two times the 
CDF between the two frequencies, or two harmonic frequencies that are separated by 
the smallest harmonic of the control signal (e.g. CDF or CDF+ SCDF. If more than 
one frequency pair is located then the frequency pair having the greatest separation 
from other harmonics is provided to the speed estimation mechanism 622. 

The speed estimation mechanism 622 estimates the speed of the motor 104 based on 
the selected candidate frequency pairs remaining after the CDF harmonics have been 
removed. The speed estimation is based on a comparison of the remaining selected 
candidate frequency pairs with a known mathematical model relating the induced 
frequencies in the motor drive connections 112 with the motor speed. If multiple 
corresponding frequency pairs have been located then the speed estimate obtained 
from each of these pairs is averaged. 

Embodiments of the present invention may be implemented in any conventional 
computer programming language. For example, embodiments may be implemented 
in a procedural programming language (e.g. "C") or an object oriented language (e.g. 
»C++"). Further embodiments of the invention may be implemented as pre- 
programmed hardware elements, other related components, or as a combination of 
hardware and software components. 
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Embodiments can be implemented as a computer program product for use with a 
computer system. Such implementation may include a series of computer instructions 
fixed either on a tangible medium, such as a computer readable medium (e.g. a 
diskette, CD-ROM, ROM, or fixed disk) or transmittable to a computer system, via a 
modem or other interface device, such as a communications adapter connected to a 
network over a medium. The medium may be either a tangible medium (e.g. optical 
or electrical communications lines) or a medium implemented with wireless 
techniques (e.g. microwave, infrared or other transmission techniques). The series of 
computer instructions embodies all or part of the functionality previously described 
herein. Those skilled in the art should appreciate that such computer instructions can 
be written in a number of programming languages for use with many computer 
architectures or operating systems. Furthermore, such instructions may be stored in 
any memory device, such as semiconductor, magnetic, optical or other memory 
devices, and maybe transmitted using any communications technology, such as 
optical, infrared, microwave, or other transmission technologies. It is expected that 
such a computer program product may be distributed as a removable medium with 
accompanying printed or electronic documentation (e.g. shrink wrapped software), 
preloaded with a computer system (e.g., on system ROM or fixed disk), or distributed 
from a server over the network (e.g., the Internet or World Wide Web). Some 
embodiments of the invention may be implemented as a combination of both software 
(e.g. a computer program product) and hardware (termed mechanisms). Still other 
embodiments of the invention may be implemented as entirely hardware, or entirely 
software (e.g. a computer program product)." 

It is apparent to one skilled in the art that numerous modifications and departures 
from the specific embodiments described herein may be made without departing from 
the spirit and scope of the invention. 
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